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Abstract 

Meixner oscillators have a ground state and an 'energy' spectrum that is equal- 
ly spaced; they are a two-parameter family of models that satisfy a Hamiltonian 
equation with a difference operator. Meixner oscillators include as limits and par- 
' ticular cases the Charlier, Kravchuk and Hermite (common quantum-mechanical) 

harmonic oscillators. By the Sommerfeld- Watson transformation they are also re- 



<^ , lated with a relativistic model of the linear harmonic oscillator, built in terms of 

' the Meixner-Pollaczek polynomials, and their continuous weight function. We con- 

■ struct explicitly the corresponding coherent states with the dynamical symmetry 

, group Sp(2,3f?). The reproducing kernel for the wavefunctions of these models is 

' also found. 



PACS: 03.65. Bz, 03.65.Fd 



1 Introduction 



^ . An oscillator is called harmonic when its oscillation period is independent of its energy. 

In quantum theory, this statement leads to its characterization by a Hamiltonian operator 
whose energy spectrum is discrete, lower-bound, and equally spaced 

Hipn = Eni'n, -En ~ ^ + coustaut, n = 0, 1, 2, . . . . (1) 

Within the framework of Lie theory, this further indicates that only a few choices of 
operators and Hilbert spaces are available if the Hamiltonian operator is incorporated 
into some Lie algebra of low dimension. 

If we relax the strict Schrodinger quantization rule, we find a family of harmonic 
oscillator models characterized by Hamiltonians that are difference operators (rather than 
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differential operators). Their spectrum is also (|1|), with n either unbounded, or with 
an upper bound A^. The wavefunctions are still continuously defined on intervals either 
unbounded or bounded, but the governing equation will relate their values only at discrete, 
equidistant points of space; the Hilbert space of wavefunctions will then also have discrete 
measure. Thus 'space' appears to be discrete rather than continuous. The two-parameter 
family of Meixner oscillator models, to be examined here, is harmonic. Limit and special 
cases of Meixner oscillators will be shown to include the Charlier, Kravchuk, and the 
ordinary Hermite quantum harmonic oscillator models. 

The Hermite, Charlier and Kravchuk oscillator models are reviewed in Section 2, 
together with their limit relations. Their common salient feature is to possess associated 
raising and lowering operators for the energy quantum number n in (||). In the first two, 
Hermite and Charlier, the Hamiltonian operator further factorizes into the product of 
these raising and lowering operators; the relevant Lie algebra is the Heisenberg-Weyl one, 
which can be extended to the two-dimensional dynamical symplectic algebra sp(2, 3?) = 
so(2,l) = su(l,l). 

The Meixner oscillator model is introduced in Section 3, using well-known properties 
of the Meixner polynomials and its difference equation. The dynamical symmetry is also 
sp(2, 3ft). It is then natural to build the coherent states of Meixner oscillators in Section 4. 
Section 5 establishes the analogue for Meixner wavefunctions of the well-known property of 
the Hermite functions to self-reproduce under Fourier (and fractional Fourier) transforms. 
This property is of great interest in the processing of signals by optical means . Section 
6 shows that limits and special cases of the Meixner oscillator are the Hermite, Charlier, 
and Kravchuk oscillators. In this Section we discuss also how the Meixner oscillator is 
related the radial part of the nonrelativistic Coulomb system in quantum mechanics and 
a relativistic model of the linear harmonic oscillator, built in terms of the continuous 
Meixner-PoUaczek polynomials. 

2 Hermite, Charlier and Kravchuk oscillators 

In this Section we collect for the reader the basic facts on the Hermite, Charlier and 
Kravchuk oscillators. We introduce their Hamiltonian operator and wavefunctions, as 
well as raising and lowering operators for each oscillator model. Finally, we show the 
limit relations whereby Charlier and Kravchuk 'discrete' oscillators converge both to the 
quantum-mechanical (Hermite) harmonic oscillator. 

2.1 Hermite (quantum- mechanical) oscillator 

The linear harmonic oscillator in nonrelativistic quantum mechanics is governed by the 
well-known Hamiltonian 

H'^iO = ^ie~ dl) = hu [a+a + 1/2] , (2) 
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where ^ = /hx is a dimensionless coordinate (m is the mass and uj is the angular fre- 
quency of a classical oscillator); we indicate = d/dC,, and the creation and annihilation 
operators are defined as usual: 

a+ = -^(e-9^), a=^{^ + d^), K a+] = l. (3) 

Eigenfunctions of the Hamiltonian (0) are expressed in terms of the Hermite polyno- 
mials Hn{C,), n = 0,1,2, ... . Their explicit form is 

HniO = 2Fo{-\n, \[l - n]; -1/^^) = n\ ^ \,[^ y^^^y^ , (4) 

where [n/2] is or ^{n — 1) according to whether n is even or odd. Hermite polynomials 
are orthogonal and of square norm c„ under integration over ^ G 3?, with measure dC,, 
where 

p"(0 = e~^', c„ = v^2"n!. (5) 



Therefore, the normalized wavefunctions 
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C(0 = Vp"(0/Cn H^niO = I ^n(0 e~^", 71 = 0,1,2,... , (6) 



are orthonormal and complete in the Hilbert space £?{^), commonly used in quantum 
mechanics, namely 

/oo 
dim) m) = Sn,k , E c(o m') = sic - a ■ (t) 

Their corresponding eigenvalues under (|^) define the energy spectrum of the harmonic 
oscillator, and are ([^) in the form 

En = huj{n + l). (8) 



2.2 Charlier oscillator 

A difference (or discrete) analogue of the linear harmonic oscillator (|]), can be built on 
the half-line in terms of the Charlier polynomials Cn{x;n), for any fixed /i > and 
n = 0, 1, 2, . . . 0. Charlier polynomials are defined as |1 

Cn{x; /.) = 2 Foi-n, -x; /i-^) = ± (9) 

where (a)„ = r(a + n)/r(a) = a{a + 1) ■ ■ ■ (a + n — 1) is the shifted factorial and T{z) is 
the Gamma function. 
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The Hamiltonian for the Charher oscillator model is a difference operator | 
where by definition 



(10) 



eiy9^f^a:) = f{x±y) (11) 

is a shift operator by y with the step /ii = l/y/2jl . Eigenf unctions of ([10|) have the same 
eigenvalues (§); they are orthogonal with respect to the weight function 

and have the form 



no = ^ + ^/^i) Cn(/i + /i) • (13) 

It is clear from the definition (|^) that the Charlier polynomials are self-dual: C„(x; /i) = 
Cx{n; /i); therefore the Charlier functions (0) satisfy two discrete orthogonahty relations 

oo oo 

E ^^i^k) n^k) = 5m,n , E V'."(en^) ^."(^n) = 5n,m, (M) 

fc=0 k=0 

where = (k — n)hi. These are the discrete analogues of the continuous orthogonality 
and the completeness relations in (|^). 

As in the nonrelativistic case (Q), it is possible to factorize the Hamiltonian ([T0| ) 

H''{i) = ?ioj{h+h+\), (15) 

by means of the difference operators 



6 = + 1 + i/h^ e^^^^ b+ = ^f, + C/hi e-'^^^^ - . (16) 

These operators satisfy the Heisenberg commutation relation 

[6,6+] = 1, (17) 
and their action on the wavefunctions (p!3| ) is 
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2.3 Kravchuk oscillators 

Another difference analogue of the harmonic oscillator |^ can be built on the finite interval 
[0, A^], where is some positive integer, in terms of the Kravchuk polynomials § 

K^ix;p,N) = ,F,i-n, -x; -N-p'') = ± ^'f'^'^'l' . (19) 

This is a family of polynomials, parametrized by < p < 1, of degree n = 0, 1, 2, ... in 
the variable x. 

The corresponding Kravchuk oscillator Hamiltonian is a difference operator with step 

h2 = v/2Ar^ §, 

H^{i) = hu:\2p{l-p)N + \ + {\-p)i/h2 

- ^p{l-p) [a(0 e'^^^^ + a{i - h) e'^^^^] | , (20) 

a(0 = ^JiqN-^/h2){pN+l + ^/h2). 

The energy spectrum is the same as in (^, except that in the Kravchuk case there are 
only a finite number of energy levels n = 0, 1, A^. 

The Kravchuk polynomials are orthogonal with respect to the binomial measure 

pK(x)=0"g^-^ C% = N\/r{x + l)r{N-x + l). (21) 

The eigenfunctions of the difference operator (|20|) are 



C(0 = {-ir^C-N{j^ypHpN + ^/h2)K4pN + ^/h2;p,N), (22) 

where is the binomial coefficient. The Kravchuk polynomials (|19]) are also self-dual, 
and therefore the Kravchuk functions (|2^) satisfy the discrete orthogonality and com- 
pleteness relations over the points C,j = {j — pN)h2: 

N N 

E (e,) = , E ^l^^n) ^"^{ik) = 5n,k , (23) 

i=o i=o 

for = 0, 1,...,A^. 

Now, it has been shown in H that the difference operators 



A{i) = {l-p)a{i)e^^'^-pe-^^'^a{i) + ^p{l-p)[{2p-l)N + 2i/h2], (24) 

A+(0 = (1 - p)e-'^^^« a(0 -pa{i) e'^^'^ + ^p{l - p) [{2p - l)Ar + 2^7/12] , (25) 
together with the operator 

Ao{0 = ^[H''{0-\{N + l)\ , (26) 



close under commutation as the algebra so(3) of the rotation group, 

[Ao, A] = -A, [Ao, A+] = A+, A] = 2Ao. (27) 

The action of the operators (p^-p5D on the wavefunctions ( p2D is given by 

AiOi^'^iO = ^n{N-n + l)^l_,{i), (28) 

^^(0C(0 = v/(^ + l)(iV-n)C+i(0- (29) 

We note that the Kravchuk oscillator was applied recently in finite (multimodal, shallow) 
waveguide optics 0. 

2.4 Limiting cases 

Among the previous models, the Kravchuk oscillator is the most general; it limits to the 
Charlier oscillator; in turn, the latter limits to the common Hermite harmonic oscillator 

i: 

Kravchuk — > Charlier. Because of the limit relation |^ 

lim Kn{x- fx/N, N) = Cn{x; /i) (30) 

A"— ►CXD 

between the Kravchuk ( |T9|) and Charlier (P) polynomials, when N —>■ oo and p = fi/N —>■ 
0, the operators H^{^), A{S,)/\/N and A'^{^)/\/N reduce to the Charlier Hamiltonian 



(p!5|), and the lowering and raising operators ([iq) for the Charlier functions, respectively. 
The so(3) algebra ( P?] ) in turn contracts to the Heisenberg-Weyl algebra ([T^). 
Charlier — > Hermite. In the limit when the Charlier parameter /i tends to infinity, 
we have 

lim /ir" C„(/i + i/h^- /i) = (-1)" Hn{i) . (31) 



Similarly, in the limit /i — » oo, the operators (|T6]) become b ^ a, a^, and H'^{^) — >■ 

H^{^). The Charlier functions (0) coincide then with the Hermite functions (^), i.e. 

lim h^'^' r^iO = mo ■ (32) 
Kravchuk — > Hermite. From the limit relations H 



hm (-1)" ^Jc%ip/qr K^aipN + ^/h,; p, N) = -== iJ„(0 , (33) 

and 

hm h,' p^pN + i/h^) = e-^' , (34) 



it follows that 

hm h,'/'mO = MO- (35) 

Also, when oo, the operators H^{^), A{^)/\/N and A~^{^)/\/N reduce to the 

Hermite Hamiltonian (^, annihilation a(^), and creation a~^{C,) operators (0) for the 
ordinary quantum harmonic oscillator, respectively. The so (3) algebra (^) of this finite 
oscillator contracts to the Heisenberg-Weyl algebra of quantum mechanics. 
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3 Meixner oscillators 

We now organize the properties of the Meixner polynomials ^ according to the scheme 
followed in the previous Section. Known orthogonality relations for the Meixner polynomi- 
als lead to orthonormal functions and a difference Hamiltonian operator, whose spectrum 
is the set of energy levels (|lD. 

3.1 Meixner polynomials and functions 

The Meixner polynomials Q are Gauss hypergeometric polynomials 

M„(e; /5, 7) = 2Fi(-n, -e; /3; 1 - I/7) = Mi:{n; /5, 7) • (36) 

They form a two-parameter family of polynomials, for j3 > and < 7 < 1, of degree 
77, = 0, 1, 2, . . .. Their orthogonality relation is 

00 

^ p"(m)M„(m;/5,7)Mfc(m;/5,7) =rf„5nfc (37) 

m=0 

with respect to the weight function and square norm 



Hence, the wavefunctions of the form 

C(e; A 7) = (-l)Vp"(ORM„(e; /3, 7) , (39) 

satisfy the discrete orthogonality relations 

00 00 

Y: C(e; P, 7) C(e; P, 7) = Sn,k , E C(e; P, 7) C(r; 7) = ^c,^ , (40) 

5=0 n=0 



as a consequence of (|37D and the self-duality of Meixner polynomials (|36|). Henceforth 
we shall supress for brevity the super-index M from all operators and functions of the 
Meixner oscillator model. 

The Meixner polynomials ( p6D satisfy the three-term recurrence relation |^ 

[n + (n + /?)7 - (1 - 7)e]M„(e; P, 7) = (ri + PhM,,+,{C; P, l) + nM^_,{^- P, 7), (41) 

and the difference equation in the real argument 

[7(e + P)e^' + ^e-^' - (1 + 7)(e + P/l) + {l + l){n + /3/2)] M^{^- /?, 7) = . (42) 
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Hence, the functions ( |39| ) are eigenfunctions of the difference Meixner Hamiltonian oper- 
ator 



H{0 



1 — 7 ^ 1 — 7 



/i(Oe^« + - l)e 



with eigenvalues 



n = 0,l,2,... . 



(43) 
(44) 

(45) 



En = n+^(3, 

3.2 Dynamical symmetry algebra Sp(2,§f^) 

As in all previous cases, we can construct the dynamical symmetry algebra (see, for 
example, @, |^) by factorizing |ll[] the difference Hamiltonian (|^). Indeed, one can 
verify that 

H{0 = BB+ + ll3-l, (46) 
where B = B[C,) and B~^ = B^{C,) are the difference operators 

1 



B 



B- 



1 



e + le^^«- V7(e + /5-l)e 



(47) 
(48) 



a/I - 7 V 

It is essential to note that the factorization of the Hamiltonian (|^), in contrast to the 
case of the harmonic oscillator (H) and the difference model (|T5|) , does not lead immediately 
to a closed algebra consisting of if(^), B and B~^. To obtain such an algebra, we compute 
the explicit form of the commutator between the last two. 



[B,B' 



+ ^(e^^^/^4^M0e^^^ + e-^^e^(0 J^^^e-i^O ■ (49) 



1-7 V" V e + "^^^V ^ + P 

The right-hand side of (^9]) suggests that it is necessary to introduce new operators 



C 



1 




7 


7 


1 




7 


7 



77 



e + le^^«5+ = ^-/i(ne^^ 
\/7 



(50) 



(51) 



These new operators have the following commutation relations with the Hamiltonian 
operator 



-C + ^iH+l-l(3), 
V7 



C+- (H + l-lP). 



(52) 
(53) 
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We now build the difference operators 

7 



(54) 



— Ki)e'' + 7^ e-^^/.(0 - + 1/3) . (55) 



1 — 7 1 — 7 1 — 7 

= C-^(iJ+l-i/5) 
V7 



1 — 7 1 — 7 1 — 7 

Together with Kq = H, they now form the closed Lie algebra sp(2,3?), 

[Ko,Ki] = ±Ki, [K.,K+] = 2Ko. (56) 

The raising and lowering operators and are connected by with the cartesian 
generators 

K, = -zl{K^-K^) = zi[;.(Oe«e-e-V(0], (57) 
K2 = -UK+ + K^) = -^^^[M0e^^ + e-S(0] + ^(e + |/5). (58) 

The invariant Casimir operator in this case is 

= Kl - Kl - Kl = Kl -Ko- K+K^ = i/5(|/? - 1)/ . (59) 

The eigenvalue — 1/3 of the Casimir operator determines that the model realizes the 
unitary irreducible representation D~^{—P/2) of the Sp(2,3?) group. The eigenvalues of 
the compact generator Kq{^) in such representations are bounded from below and equal 
to lf3 + n, n = 0, 1, 2, . . .. In other words, a purely algebraic approach enables one to find 
the correct spectrum of the Hamiltonian H{^) = -ft'o(0 (i^). 

The action of the raising and lowering difference operators and on the wave- 
functions (|39D is given by 



^+^n(^;/3,7) = '«n+l^n+l(^;/5,7) , K_llJn{^;f3,-f) = KniJn-l{^;P,l), (60) 



where Kn = yn{n + (3 — 1). Hence the functions ?/'„(^;/5,7) can be obtained by ra-fold 
application of the operator to the ground state wavefunction 

= -r^—KiM^;P,7), (61) 



^o(e;A7) = V(1-7Fp(0- (62) 
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3.3 Unitary equivalence in the second parameter 

Observe that the eigenvalues of the Casimir operator as well as the matrix elements 
(|60|) of the operators and K_, do not depend on the second parameter, 7, of the 
Meixner wavefunctions. Therefore the basis functions corresponding to two distinct 
values of the parameter 7, must be intertwined by a unitary transformation. To find its 
explicit form we may compare two sets of the generators Kq^ Ki and K2 [see formulas 
(^) and (|57| , |58D ], corresponding to different values 7 and 7'. 

Introducing angles 9 and 9' such that 7 = tanh^ 7' = tanh^ and 6 = 9' — 9, the 
relation between the two sets of generators is written as 

Kq = cosh 6 K'q + sinh 6 K'2 , 

K, = K[, (63) 
K2 = sinh 6Kq + cosh 5 . 

This shows they are related by a boost in the — 2 plane by the hyperbolic angle 6 G 
3?. Consequently, the wavefunctions (pOf ) with different values of the parameter 7 are 
connected by 

00 

MC, P, 7) = e'""'^' M^; P, 7) = E KCk MC, P, 7) • (64) 

fc=0 

The last expression is the matrix form, with elements 

00 

5=0 

= (tanh i5)"+'= (cosh IS)-^ Mn{k, /?, tanh^ ^S) . (65) 

In deriving ( |65D we have used the addition formula for the Meixner polynomials (^) given 
in Ref. Eq. (A. 6). 



4 Coherent states 



The dynamical symmetry of the Meixner oscillator model (P3|), allows us to construct 
two kinds of coherent states [13, 14]. Recall that in the case of harmonic oscillator (|^) 
coherent states are defined as eigenstates of the annihilation operator a(^) []T3[. Coherent 
states for the model (|43D can be defined either as eigenstates [^ of the lowering operator 
K_{^), or by acting on the ground state (|62D with the operator exp[^_R'+(^)] [Q. This 
gives rise to two distinct coherent states. 
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4.1 The Barut-Girardello coherent states 

Characterizing the Barut-Girardello coherent states by the complex number z E C, which 
is the eigenvalue under the lowering operator, 

K^U^;P,j)=zM^;P,j), (66) 

these coherent states can be expanded in terms of the wavefunctions (p9D, 



(3,7) = 1: ^— M^; 7) • (67) 

n=0 Jn\{(3)n 



Using the generating function p for the Meixner polynomials ( p6| 



00 j.n 



Y: ^M„(e;/?,7) = eSFi(-e;/?; -^t), (68) 



their explicit form is found to be 



,(e;/5,7) = e-^v^iFi 



7 - 1 ■ 

V7 . 



^o(e;/5,7)- (69) 



These coherent states are overcomplete and therefore nonorthogonal, 



E0:(e;/3,7)0.'(e;A7) = {z*zT-^y'T{p)i^^,{V^), m 

where Iu{z) is the modified Bessel function. 
4.2 Perelomov coherent states 

The second definition of generalized coherent states is due to Perelomov [l^; it is built 
through the action of the group operator exp(^ii"+) on the ground state ipoi^^; 

xc(e;/3,7) = (i-|Ciy/'exp(Ci^+)^o(e;A7) 



OO 



= (1-lcrr E\/^c^.(e;/3,7), (71) 

where C is a complex number such that \(\ < 1. 

Using the generating function for the Meixner polynomials 

E ^ ^" Mn{^-^ A 7) = (l - (1 - tr'~', (72) 

we find 

xc(e; A 7) = (1 - \C\y^' (^^^y^^^ V^cr^-' M^; 7) ■ (73) 

These coherent states satisfy the relation [cf. ([70|)] 

OO 

E xc(e; P, 7) xc'(e; A 7) = [(1 - ICf) (1 - (1 - CCT' ■ (74) 
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5 Reproducing transforms 



Consider the task to find a reproducing kernel for tlie Meixner functions (p9D, defined by 
tlie relation 

oo 

^ /Ci(e, Mi'. 7) = Mi: 7) • (75) 
The quantum mechanical analogue of this expression is the property of Hermite functions 



to reproduce under fractional Fourier transforms of angle r for t 
Fourier transform of kernel exp(i,^,^') corresponds to r = ^Tr 



e*"^; the common 
T7| . The finite-difference 
Fourier- Kravchuk transform has the same property on the Kravchuk functions, and has 
been shown recently to apply to shallow multimodal waveguides with a finite number of 
sensors 0. 

Using the dual orthogonality relation of the Meixner functions (|40|) , the explicit form 
of the kernel JCt{i,i') is found for \t\ < 1, 



(76) 



n=0 



It is a bilinear generating function for the Meixner functions. By the definition (|76D , 
the reproducing kernel JCtid') is symmetric with respect to exchange of ^ and and 
because of the orthogonality relation (^) it has the property 



(77) 



Reproducing kernels for the Charlier (|T2D and Kravchuk (|2^) functions have been discussed 
in ITH], whereas the cases of the g- Hermite and Askey- Wilson polynomials have been 



considered in and |2y, ^ , respectively. 
Substituting (|39|) in (^), we can write 



/c.(e,n = vp(0p(e)(i-7)^E 



n=0 



(7t)"M„(e;/3,7)M.(r;/3,7) 



(78) 



The sum over n in ( |T5D is the bilinear generating function (Poisson kernel) for the Meixner 
polynomials [Q, 



^ iq!^rM„(e;A7)M„(e';A7) 



n=0 



= (l-t)-/5-«-«'(l-t/7)«+«'2Fi 
Thus the kernel JCt{i,i') is written as 



,-;^^ (l-7)^(l-t)^+^' p 



1 _ ^t)^H'+P 



^(1-7^) 
(t-7)2 

t(l-7^) 



7(1 -t)2. 



(79) 



(80) 



12 



For integer ^ and ^' we have the hmit 

£m/C,(e,0=%?'- (81) 

In this hmit, the relation (^) coincides with the dual orthogonality (^Op of the Meixner 
functions. 

The limit of (80) when t ^ i {t ^ ^n) corresponds to a discrete analogue of the 
classical Fourier-Bessel transform; whereas the latter integrates over the nonnegative half- 
axis, the former sums over the integer points ^ = 0, 1, . . .. The limit is 



lim/Ci(e,0 



/c.(e,0 

(1 - i7)?+?'+/3 

It is easy to verify that for integer ^ and 



— — - (-2.)«+^')/^(l-7)^ 
P(0 ) T-T7T77T^^ 2F1 



7)^ 



27 



f82) 



5'=o 



(83) 



6 Limit and special cases 

The difference model of the Meixner linear harmonic oscillator family ( ^Bj ) contains as 
limit and particular cases all the models of Hermite (0), Charlier (0) and Kravchuk (pO]). 
We make these limits explicit below. Here we discuss also the corresponding relations 
with the radial part of the nonrelativistic Coulomb system in quantum mechanics and 
a relativistic model of the linear oscillator, built in terms of the continuous Meixner- 
PoUaczek polynomials. 



6.1 Meixner — > Hermite 

From the recurrence relation for the Meixner polynomials (0), it can be show that the 
following limit to the Hermite polynomials holds: 

hm (2z/)"/2 M„ ( ^ + ^^ . ^\ = (_i)n (84) 
f — »oo \ — 7 7 / 

Furthermore, measures and normalization coefficients relate as 

limV2i;(l-,)''-V(^^) = (85) 

lim (2z/)"/2 (1 _ ^)'^/27^^ = 72^. (86) 
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The wavefunctions ( pQ] ) with argument (z/ + V2z/,^)/(l — 7) and (5 = vj^^ coincide in 
the hmit v 00 with the wavefunctions of the hnear harmonic oscillator (1^), z.e., 



lim 



; -,7 



€(0 • 



7) 



The combination -ft'o(0 — ^/'^l reproduces, in the same limit, the product a"*" (0*^(0 



whereas the matrix elements of ^J'~) /vK±{^) converge to the creation and annihilation 
operators a^{C) and a(^), respectively. The Meixner oscillator family (|43|) thus contains 
as a limit case the linear harmonic oscillator (ET) of quantum mechanics. 



6.2 Meixner — > Charlier 

It is known that the Meixner (|36D and Charlier @ polynomials are connected by the limit 
relation |, § 

hm M„(e;/3,/i//3) = C„(e,/i). (88) 

Hence in the limit when (3 ^ 00 and 7 = — > 0, from (^) one obtains the reproducing 
kernel for the Charlier functions HTHi: 



/c^(e,n 



hm m.o = e 



-2i/i) 



2-ro 



2/i 



Using the limit relation (p8| ) it is easy to check that 

lim V'„(/i + e//ii;A/i//3) = V'n(0 

p— >oo 



^9) 



(90) 



Hence the wavefunctions ( |55| ) with argument + C/^i and parameter 7 = coincide, 
in the limit when /3 — >^ 00, with the wave funcitons of the difference (discrete) model of 
the Charlier oscillator (0). In the same limit, the combination H{s) + |(1 — /3), where 
s = jJL + ^/hi, reproduces H'^{^)/htu, whereas P~^^'^K±{^) tend to the raising and lowering 
operators and B, respectively. 



6.3 Meixner 



Kravchuk 



The Kravchuk polynomials ( |T9| ) are also a particular case of the Meixner polynomials 
(^), with the parameters /3 = — and 7 = —p/ (1 — p), that is, 

ir„(e;p, N) = M„(e; -N,p/{p - 1)) . (91) 

In this case, from (^) one obtains the reproducing kernel for the Kravchuk functions ( [2^ ) 
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]c^{^, a = \/(-2w)?+«' cic^; (1 - (1 - ^)p 



1 

2p(l-p), 



(92) 
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As follows from the relation (pT|) for f3 = —N, 7 = — — p) and argument 



pN + ^J2p{l — p)N^, the model (^) coincides with the difference model of the Kravchuk 
oscillator (EOl). 



6.4 Meixner — > Laguerre 

The limit relation 

limM„(x//i;/5,l-/i) = -—Lt\x), (93) 

where L'^{x) are the Laguerre polynomials, enables us to consider the nonrelativistic Coul- 
omb system as another limit case of the difference model (^31). Indeed, from (|39|) and (pB]), 
it follows that 



lim -^iprr{2kr/h; 2l + 2,l-h) = V2k^Rn ii2kr), (94) 



RnA^) = (-ir^ (^^^''^"_^^), ^'e-/^^^'+^(x) (95) 



where 



is the radial wavefunction of the Coulomb system (see, for example, |£4|)' radial 
variable, / and n + / + 1 are orbital and principal quantum numbers respectively, and 
k = me^/fi^in + 1 + 1). 

In the same limit h 0, the generators -^'0(2;) and K±{x), where x = 2kr/h re- 
produce the well-known generators of the dynamical symmetry algebra su(l,l) for the 
nonrelativistic Coulomb model iSd, 



rdl + a - "^—^ - kh- 



(96) 



r 

J± = -jQ + kr + {rdr + \). (97) 

This connection between the Meixner polynomials and radial wavefunctions for the 
nonrelativistic Coulomb system can be used for constructing a g-analogue of the Coulomb 
wavefunctions in terms of the g- Meixner polynomials (see Refs. P3| , pUf). 

6.5 Meixner-Pollaczek (relativistic) oscillators 

There is the family of Meixner-Pollaczek polynomials 

pVx; 0) = e-'"<^ F(-n, A - ix; 2A; 1 - e^^^^), (98) 

which satisfy the orthogonality relation 

r Pni^; <P) Pn'i^; 0) P'iO = Sn,n' (99) 
J — oo To: 
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with respect to a continuous measure with the weight 

p^iO = :^(2sin0)2^|r(A + ^Ol'exp[(20 - tt)^]. (100) 
Zn 

The reason why we mention these polynomials here is the following. In Ref. it was 
shown that the Meixner polynomials M„(^; /3,7) and the Meixner-PoUaczek polynomials 
(p^) are in fact interrelated by 

„—in(p 

(e; 0) = ^(2A)„M„(ze - A; 2A, e"^^^). (101) 
nl 

The transition from the discrete orthogonality ( P7| ) for the Meixner polynomials M„(^; /?, 7) 
to the continuous one (p9D is analogous to the well-known Sommerfeld- Watson transfor- 
mation in optics and quantum theory of scattering. 



In the relativistic model of the linear harmonic oscillator, proposed in the wave- 
functions in configuration space are expressed in terms of the Meixner-Pollaczek poly- 
nomials ( pHD and their weight function ( p.00| ) with the specific value of the parameter 
(f) = ^TT. The same model in the homogeneous external field gx corresponds to the value 
of the parameter (p given by arccos((7/mca;), where m and u have the same meaning as 



in the classical case, and c is the velocity of light [pSj |29|| . In other words, the relation 
( |101| ) gives the connection between the relativistic harmonic oscillator and the Meixner 
oscillator, discussed in this paper. 
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